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Abstract: We perform numerical simulations of the gravitational collapse of a k-essence 
scalar field. When the field is sufficiently strongly gravitating, a black hole forms. However, 
the black hole has two horizons: a light horizon (the ordinary black hole horizon) and a 
sound horizon that traps k-essence. In certain cases the k-essence signals can travel faster 
than light and the sound horizon is inside the light horizon. Under those circumstances, 
k-essence signals can escape from the black hole. Eventually, the two horizons merge and 
the k-essence signals can no longer escape. 
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1 Introduction 

It has been experimentally verified that the universe has recently entered an era of accel- 
erated expansion. This discovery has led physicists to pose the coincidence problem: why 
does the accelerated expansion of the universe coincide with the onset of matter domina- 
tion? k-essence [1, 2] has been proposed in order to answer this question, k-essence models 
postulate the existence of a scalar field with non-standard kinetic terms that allows for the 
field to have negative pressure once the matter dominated era begins. 

It has been shown that if k-essence is to solve the coincidence problem, then fluc- 
tuations of the k-essence field must necessarily propagate superluminally at some point 
[3]. However, it has been argued that as long as causality is not violated this should be 
physically acceptable [4]. 

Given the possibility of superluminal propagation, it is interesting to ask how a k- 
essence field would act during gravitational collapse. It has been found that the stationary 
solution describing k-essence scalar field accretion into a background Schwarzschild black 
hole allows for the existence of two horizons, one corresponding to the fluctuations of the 
k-essence scalar field and the other being the normal light horizon [5] . This is interesting as 
the existence of two horizons also occurs in other theories that allow superluminal motion 
[6, 7]. It is possible that the stationary solution allowing for the existence of two horizons 
may be avoided kinematically given reasonable initial conditions. It is also possible that 
once back reactions of the scalar field on the metric are considered the two horizons may 
reduce to one. Another possibility is that gravitational collapse results not in the formation 
of a black hole but a naked singularity. 

In order to see which of the above mentioned possibilities are realized during the 
gravitational collapse of a k-essence scalar field, one can perform numerical simulations of 
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this process. Such simulations have often been done before for an ordinary scalar field [8, 9] 
and the methods used can be modified and generalized to treat k-essence. In this particular 
study, we will numerically simulate the self-consistent gravitational collapse of spherically 
symmetric scalar fields described by several different k-essence Lagrangians. The knowledge 
of what k-essence models allow for what type of horizons may help restrict the parameter 
space of acceptable k-essence models. Numerical simulations may also indicate whether or 
not there are any unique experimental signatures associated with the creation of a black 
hole from a k-essence scalar field. 

The equations and numerical methods used are described in section 2 and the appendix. 
Results are presented in section 3 and conclusions discussed in section 4. 



2 Methods 

The action for self gravitating k-essence is 



/ ^ d ' X + (2 - 1} 



Here is the spacetime scalar curvature, k = 8irG where G is Newton's gravitational 
constant, and X = — ^V a </>V a </> where (j> is the k-essence scalar field. Varying the action 
with respect to the metric yields the Einstein field equation 

G a b = K T a b (2.2) 

where G a b is the Einstein tensor, and the k-essence stress-energy tensor T a b is given by 

Tab = £xV a <AV fe + Cg ah (2.3) 

Here, we introduce the shorthand Cx to mean dC/dX (and later we will also use Cxx to 
mean d 2 C/dX 2 ). Varying the action with respect to the k-essence scalar field 4> yields the 
k-essence equation of motion 

r b V a V^ = (2.4) 

Here g ab is an effective inverse metric associated with the k-essence scalar field and is given 
by 

g ab = C x g ab - £xxV a <f>V b <f> (2.5) 

We perform numerical simulations of equations (2.2) and (2.4) in the case of spherical 
symmetry. For an ordinary scalar field, such simulations are often performed using polar- 
radial coordinates [8]. Here the radial coordinate is chosen to be the usual area radius, 
while the time coordinate is chosen to be orthogonal to the radial coordinate. However, 
polar-radial coordinates only last until the first trapped surface forms. This would not be 
suitable for a simulation of k-essence gravitational collapse because k-essence travels faster 
than light and thus is not trapped by ordinary trapped surfaces but instead by the trapped 
surfaces of the effective metric of eq. (2.5). Thus a k-essence black hole has not yet formed 
when the first trapped surface forms and a simulation using polar-radial coordinates stops. 
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Figure 1. <f> vs r at the initial time i = 0for J 4 = 0.15 



Instead, we use a different method with a set of coordinates that includes a region within 
the black hole horizon. This method uses maximal slicing, and a radial coordinate based 
on length rather than area, and is described in detail in the appendix. 

Another complication of k-essence collapse has to do with the surfaces of constant 
time. In an evolution of the Einstein field equations, the surfaces of constant time must be 
spacelike. That is, we must have g u < 0. For an ordinary scalar field, this condition is also 
enough to render the evolution of the scalar field well posed. However, k-essence evolves in 
an effective metric whose light cones are different from those of the spacetime metric. In 
order that the surfaces of constant time are also Cauchy surfaces for the k-essence field, we 
must have g < 0. In the simulations, it is sometimes the case that g — > without any 
singularity forming and sometimes even without the formation of any trapped surface. In 
this case we are forced to stop the simulation at that point without obtaining a clear answer 
as to the outcome of the collapse process. This also suggests the disturbing possibility that 
the light cones of the k-essence field could tip so much that no surface could be timelike 
for both the metric and the k-essence field. In that case, the theory itself would cease to 
be well defined. 

In any case, our results are obtained by choosing a variety of k-essence Lagrangians 
C and a variety of initial data, and then evolving with the maximal slicing method, and 
noting the outome of the collapse process. 



3 Results 

We now present the results of the collapse simulations. We chose the scalar field to initially 
have dt4> = and the form: 

4> = -2A r -^eM^^) (3-D 

where A, a, and tq are constants. For all results and figures shown we used the values 
a = 1.5, and vq = 25 (see fig 1 for which A = 0.15). 
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3.1 Sound horizon Lagrangian 

We first consider a k-essence model proposed in [5] with a Lagrangian density given by: 



We chose A = for all of our simulations. In [5] it was found that when there is a steady 
flux of scalar field described by this Lagrangian falling into a Schwarszschild black hole, 
there exists a stationary solution which can allow scalar field signals to escape from the 
light horizon. This is because the k-essence scalar field can propagate at speeds different 
from that of the speed of light. Thus there is a sound horizon associated with the scalar 
field which can differ from the traditional light horizon. When this sound horizon is located 
within the light horizon, signals can escape the black hole. 

Before presenting the results, it is helpful to recall how this sort of initial data would 
evolve for an ordinary (non k-essence) massless free scalar field, both without and with 
gravity. For an ordinary non-gravitating scalar field, initial data of the type given in eq. 
(3.1) separates into an outgoing pulse and an ingoing pulse. The ingoing pulse undergoes 
interference in a region near r = and becomes an outgoing pulse, so that at sufficiently 
late times the scalar field consists of two outgoing pulses. For an ordinary gravitating 
scalar field, when the amplitude A is sufficiently small, the self-gravity of the field is weak, 
and the behavior is essentially that of a non-gravitating scalar field. However, when A 
is sufficiently large, as the ingoing pulse nears the center, it becomes trapped by its own 
self-gravity and forms a black hole. 

For the model of eq. (3.2) when a 2 is large compared to X, the Lagrangian approaches 
that of a massless free scalar field. Therefore we would expect that for large enough a we 
would get the same sort of behavior as for an ordinary (non k-essence) scalar field. To 
see whether this expectation is right, we performed simulations with a = 1.0 and with 
A = 0.05 and A = 0.15. The result of the simulation with A = 0.05 is shown in figure 
(2). Here the scalar field <f> is graphed as a function of r at time t = 45. In this amount of 
time the initially outgoing pulse has progressed to steadily larger r with approximately 1 jr 
falloff in amplitude, while the initially ingoing pulse has undergone interference near the 
center and become an outgoing pulse which has moved away from the center undergoing 
its own approximately 1/r falloff in amplitude. 

In contrast, figure (3) shows the result of the simulation with A = 0.15 (and a = 1.0). 
Again the scalar field (j) is graphed as a function of r at time t = 45. Once again, the 
outgoing pulse has progressed to steadily larger r with a corresponding falloff in amplitude. 
However, now the ingoing pulse has not become outgoing and instead remains trapped 
near the center by its own self-gravity. Indeed we would expect the ingoing pulse to not 
merely "remain trapped" but rather to undergo complete gravitational collapse and form a 
singularity at the center of a black hole. The fact that it has not done so in this simulation 
is an artifact (and one of the main features) of the maximal slicing method. In maximal 
slicing, the lapse function N becomes small in regions that are strongly gravitating. This 
slows the evolution down in those regions and prevents the simulation from reaching the 
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singularity (which would force the simulation to halt). Indeed such a "collapse of the lapse" 
is usually taken as an indication of black hole formation. This phenomenon is illustrated 
in figure (4) where the lapse for this simulation is plotted. 

Nonetheless, we would like a more precise indication of black hole formation. In general 
relativity, such an indication is provided by the formation of a marginally outer trapped 
surface. Here a two dimensional surface has an outgoing null geodesic normal vector l a and 
the surface is marginally outer trapped if V a £ a = 0. This condition simplifies considerably 
in the case of spherical symmetry, where it becomes equivalent to £ a V a R = where R is 
the area radius (i.e. A-kR 2 is the area of the spheres of symmetry). That is, we want to 
know whether there is an outgoing vector l a such that g a b£ a £ b = and £ a V a R = 0. For 
the coordinate system we use (eq. (A.l)) we find after some straight forward algebra that 
this condition becomes: 

-RK; + R' = (3.3) 

Thus, to check for the presence of a marginally outer trapped surface (and thus the for- 
mation of a black hole) we only need to compute the left hand side of eq. (3.3) and see 
whether it vanishes somewhere. 

The condition of eq. (3.3) works for the gravitational collapse of an ordinary scalar 
field. However, since k-essence can travel faster than light, we want to know not only 
whether an object forms that can trap light, but also whether an object forms that can 
trap k-essence. Since k-essence satisfies a wave equation using the effective inverse metric 
g ab , it follows that the direction of propagation of k-essence is a vector that is null with 
respect to the effective metric which satisfies g~ b l g bc = S c a - Straightforward calculation 
using eq. (2.5) then shows that 

9~b = 7-9ab + <2^prV a <l>V b <l> (3.4) 
>->X J~x 

where 



(3.5) 



C x + 2XC XX 

is the speed of propagation of the k-essence waves for this model, (we will also refer to c s 
as the sound speed). Thus the condition for a sound horizon is just that for a light horizon, 
but with the effective metric replacing the metric. That is, the existence of an outgoing 
vector £ a satisfying g^£ a £ b = and £ a V a R = 0. Some straightforward algebra using eq. 
(A.l) then shows that this condition is 

\rk\{c x + c 2 s C xx S 2 ) 



+R'( v / £ x (£ x -2Xc 2 £ xx ) - c 2 £ xx PS) = (3.6) 

In summary: the condition for a light horizon is the vanishing of the left hand side of eq. 
(3.3) while for a sound horizon it is the vanishing of the left hand side of eq. (3.6). 

In figure (5) are plotted the left hand sides of eq. (3.3) (light horizon, solid line) and 
eq. (3.6) (sound horizon, dotted line) at a time shortly after the sound horizon first forms 
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Figure 2. <j> vs r at time t = 45 for A = 0.05 and a = 1.0 
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Figure 3. cf> vs r at time t = 45 for A = 0.15 and a = 1.0 




Figure 4. the lapse function N vs r at time f = 45 for A = 0.15 and a = 1.0 



(t = 38.8). Note that each curve crosses zero more than once; however it is the outermost 
crossing for each curve that is the horizon. Also note that the sound horizon is inside the 
light horizon, as we would expect for k-essence that travels faster than light. A similar 
comparison done at t = 45 shows that both horizons have merged. In physical terms, the 
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Figure 5. the conditions for light horizon (solid line) and sound horizon (dotted line) at time 
t = 38.8 for A = 0.15 and a = 1.0 and C given by cq. (3.2) 



reason for this merger is that eventually all scalar field either falls into the black hole or 
escapes from it. Thus at sufficiently late times there is no scalar field on the horizon, which 
thus makes the k-essence effective metric equal to the spacetime metric in that region. 
Thus, the two horizons very quickly coincide early in the collapse process. This would 
indicate that it is only possible to send signals from inside the black hole for a very short 
time window right after a scalar wave packet falls into the black hole. In contrast, the 
reason for the persistence of two distinct horizons in the stationary solution of [5] is that 
the separation is maintained by a constant flux of scalar field through the horizon. 

Since the Lagrangian of eq. (3.2) essentially coincides with that of a free scalar field 
for large a one might expect to find behavior more uniquely characteristic of k-essence by 
performing simulations with small a. However, such simulations quickly run into difficul- 
ties. In particular, a simulation with A = 0.15 and a = 0.5 fails shortly after t = 34. At 
first sight this is somewhat puzzling, since at times shortly before the simulation fails, the 
scalar field shows no signs of the sort of rapid growth that would lead to singular behav- 
ior. However, this mystery was quickly resolved by an examination of eq. (A. 20). Note, 
that this equation involves division by the quantity Cx + CxxP 2 and an examination of 
the simulation showed that this quantity was tending towards zero at a particular spa- 
tial point just as the simulation was about to fail. One can understand this phenomenon 
geometrically as follows: from eqs. (2.5), (A.l) and (A. 10) it follows that 

~9 tt = -j^(C-x + CxxP 2 ) (3.7) 

So, it is evident that the program crashes when g u = 0. But this is precisely where the 
coordinate t fails to be timelike with respect to the k-essence effective metric, and therefore 
where the constant t surfaces fail to be Cauchy surfaces for the k-essence field equation. 
This sort of condition is familiar from general relativity, where the coordinates must be 
chosen so that g u < in order that the constant t surfaces remain good Cauchy surfaces 
for the evolution of the spacetime metric. However, when one is evolving self-gravitating 
k-essence it is necessary that the constant t surfaces be timelike with respect to both the 
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Figure 6. the conditions for light horizon (solid line) and sound horizon (dotted line) at time 
t = 38.5 for A = 0.15 and a = 0.5 and C given by eq. (3.8) 



spacetime metric g and the k-essence effective metric g . 

In order to avoid this problem, one is led to consider Lagrangians other than that of 
eq. (3.2). In particular, one could simply change the sign of coefficient of a 2 in eq. (3.2) 
to obtain 

(3.8) 

In figure (6) are plotted the results of a simulation using the Lagrangian of eq. (3.8). 
Here A = 0.15 and a = 0.5 and the plot is done at a time (t = 38.5) shortly after the 
formation of the light horizon. As in figure (5) what are plotted are the left hand sides of 
eq. (3.3) (light horizon, solid line) and eq. (3.6) (sound horizon, dotted line). As before, the 
horizons occur at the outermost point where each curve crosses zero. However, in this case 
the sound horizon is outside the light horizon, so there is no possibility of signals escaping 
from the black hole. Simulations using the Lagrangian of eq. (3.8) do not seem to have the 
problem with g u vanishing. Nonetheless, these simulations also fail at sufficiently small a. 
For example, a simulation with A = 0.15 and a = 0.1 fails at around t = 35. Here the 
problem is not g which remains negative, but rather the vanishing of c s . This is a familiar 
phenomenon from simulations of hydrodynamics. When the speed of propagation goes to 
zero, shock waves tend to form. One can still simulate such systems, but then special 
"shock capturing" numerical methods are needed rather than the simple finite differencing 
methods we use here. 

We also performed simulations of the Lagrangian given in eq. (3.2) in the absence 
of gravity. In this case, we found that for sufficiently small values of a the simulation 
still eventually fails due to g going to zero. This indicates that even the Minkowski time 
coordinate eventually ceases to be a valid global time coordinate for the emergent k-essence 
metric. As one might expect, no sound horizon was formed in this case. 

On the other hand, it is possible to create a sound horizon in the absence of gravity if 
one simulates eq. (3.8) with a sufficiently small a. As in the case where we did consider 
gravity, these simulations failed due to c s going to zero as opposed to g u going to zero. 



C(X) = a 2 
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Figure 7. the conditions for light horizon (solid line) and sound horizon (dotted line) at time 
t = 38.5 for A = 0.15 and C = 1.0 and C given by eq. (3.9) 

3.2 Lagrangians where t is a valid global time coordinate 

Due to the problems encountered using the Lagrangians of eqs. (3.2) and (3.8) we are led 
to consider Lagrangians with non-standard kinetic terms that do not allow either g a nor 
c s to vanish. One such Lagrangian is: 

CX x 

c = YTc + e ~ c T ~ l (3,9) 

where C is a constant. Like the previous Lagrangians considered, this one has the desirable 
feature that it reduces to the free scalar field Lagrangian at small values of X. Note that 
as long as C is positive, both Cx and Cxx are positive so g tt will not go to zero. Another 
requirement is that the k-essence wave equation (eq. (2.4)) be hyperbolic. The condition 
for hyperbolicity of eq. (2.4) is that Cx + 2AXxx > 0. (this condition is satisfied for the 
Lagrangians of eqs. (3.2) and (3.8)). For the Lagrangian of eq. (3.9), Cx + IXCxx has 
one global minimum with the value — j^+Uje^' Consequently, as long as we choose 
C > 0.45, the k-essence wave equation will be hyperbolic. The results of a simulation using 
this Lagrangian are plotted in figure (7). Here A = 0.15 and C = 1 and the plot is done 
at a time {t = 38.5) shortly after the formation of the light horizon. The plot shows the 
left hand sides of eq. (3.3) (light horizon, solid line) and eq. (3.6) (sound horizon, dotted 
line); and the horizons occur at the outermost point where each curve crosses zero. Here, 
the sound horizon is slightly outside the light horizon, so no signal can escape from the 
black hole. 

Another Lagrangian that doesn't result in either g or c s vanishing is 

C = X + CX 2 for X > 

C = X for X < (3.10) 

It is clear that as long as C and X are positive, g tt < and Cx + IX Cxx > 0. For 
negative X the Lagrangian is that of a free field so there are no issues in that regime. The 
results of a simulation using this Lagrangian are plotted in figure (8). Here A = 0.15 and 
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Figure 8. the conditions for light horizon (solid line) and sound horizon (dotted line) at time 
t = 38.5 for A = 0.15 and C = 1.0 and C given by cq. (3.10) 

C = 1 and the plot is done at a time (t = 38.5) shortly after the formation of the light 
horizon. The plot shows the left hand sides of eq. (3.3) (light horizon, solid line) and eq. 
(3.6) (sound horizon, dotted line); and the horizons occur at the outermost point where 
each curve crosses zero. Note that there is a discontinuity in the dotted line. This is due 
to the discontinuity at X = of the quantity Cxx- Here, the sound horizon is outside 
the light horizon, so no signal can escape from the black hole. Whether or not the sound 
horizon being outside the light horizon is a generic feature of k-essence models in which t 
remains a valid global time coordinate is yet to be answered. 

3.3 Cosmological Lagrangians 

Other models of interest are cosmological k-essence models which serve as dynamical at- 
tractors that cause the field to serve as a cosmological constant at the onset of the matter 
dominated era. In general these Lagrangians are of a factorizable form in (fi and X, that is 
they take the form: C = A{cj))B{X). Two such examples that can be found in the literature 
are [1, 2] : 

C(6, X) = 4t(-2.01 + 2^/lTX + 3 • 1(T 17 X 3 - 10" 24 X 4 ) (3.11) 

qr 

and 

£(</>, X) = ^(-2.05 + 2 v / IT7PO) (3.12) 

where 

f(X) = X - 10~ 8 X 2 + 1(T 12 X 3 - 10~ 16 X 4 

+ 10- 20 X 5 - 10~ 24 X 6 /2 6 (3.13) 

Note during our simulations with these initial conditions, X does not get much larger than 
unity, so we can neglect terms that are higher order in X with small coefficients. Thus, for 
our purposes, both these Lagrangians reduce to approximately: 

X) = ^(-1 + ^/lTx) (3.14) 



- 10 - 



The energy density of this field is proportional to As we know the universe is ap- 
proximately flat, we require that any fluctuation of this field be a perturbation on a large 
constant field. Thus, the factor of ^ in the Lagrangian will amount to a small overall 
factor. As the energy density will be negligibly small with large <f>, the gravitational effect 
will be small. Thus, both eq. (3.11) and eq. (3.12) should reduce to eq. (3.2) in flat space 
with a = 1 and a small overall rescaling factor multiplying the Lagrangian. This overall 
rescaling factor has no effect in the non-gravitational case. Simulations of eq. (3.2) in flat 
space with a = 1 result in no sound horizon but still suffer from having g tt go to zero. 

4 Conclusions 

Our results show that the gravitational collapse of k-essence scalar fields can lead to black 
holes, just as it does for ordinary scalar fields. However, the nature of these black holes is 
different. Rather than a single horizon they have two: a light horizon and a sound horizon. 
The light horizon is what we ordinarily think of as the black hole horizon. The sound 
horizon (the place where the k-essence itself becomes trapped) can either be outside the 
light horizon (if the k-essence travels slower than light) or inside the light horizon (if the 
k-essence travels faster than light). Thus, k-essence can allow signals to escape from a black 
hole. In the cases we study, the black hole is formed by the collapse of a pulse of k-essence. 
Eventually all of this pulse either falls into the black hole or escapes from it. Once that has 
happened, the two horizons have merged into one, and it is no longer possible for signals 
to escape from the black hole. One should be able to maintain the two horizons, and the 
possibility of signals escaping from the black hole, by feeding a black hole a steady stream 
of k-essence scalar field. This question is currently under study [10]. 

None of our simulations give any indication of the formation of a naked singularity. 
However, in certain cases our simulations fail not because of any singular behavior in the 
spacetime itself, but rather because our constant time surfaces fail to be Cauchy surfaces 
for the evolution of the k-essence scalar field. Since k-essence was originally proposed 
for cosmological models, one might wonder whether cosmological models might also be 
vulnerable to this sort of pathology. Here, it is helpful to remember that cosmology treats 
a situation that is homogeneous and therefore where all physical quantities depend only 
on time. For k-essence this means in particular that the quantity X is always postive 
in cosmology. In our simulations of the k-essence model of eq. (3.2) the constant time 
surfaces can fail to be Cauchy surfaces for k-essence only when X is negative. Thus, our 
gravitational collapse study brings to light certain properties of k-essence that would not be 
noticed in purely cosmological studies. We have also proposed two k-essence Lagrangians 
for which this Cauchy surface problem does not occur. It would be interesting to see if a 
Lagrangian of this sort could also be used for the usual purpose of k-essence: a natural 
cosmological explanation of dark energy. 

For those cases where the constant time surfaces fail to be Cauchy surfaces for k- 
essence, a natural solution to try would simply be to pick a different time coordinate. 
However, for evolution of self-gravitating k-essence the constant time surfaces need to be 
Cauchy surfaces for both the gravitational field and the k-essence scalar field. It is not 
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clear whether there is any time coordinate that will be able to do both jobs. If there is no 
time coordinate that can do both jobs, then that particular k-essence Lagrangian would be 
incompatible with gravity and would therefore not be a physically sensible theory. However, 
our study yields no such assertion of incompatibility because we have used only one time 
coordinate: that of maximal slicing. To resolve this issue, it is necessary to try other time 
coordinates to see whether they are compatible with both gravity and k-essence. This 
question is currently under study [11]. 

Acknowledgments 

We would like to thank Alex Vikman and Robert Mann for helpful discussions. The work 
of DG was supported by NSF grant PHY-0855532 to Oakland University. RA and RS are 
supported by a grant from the US Department of Energy. 

References 

[1] C. Armcndariz-Picon, V. Mukhanov, and P. Stcinhardt, Phys. Rev. Lett. 85, 4438 (2000) 

[2] C. Armendariz-Picon, V. Mukhanov, and P. Steinhardt, Phys. Rev. D 63, 103510 (2001) 

[3] C. Bonvin, C. Caprini, and R. Durrer, Phys. Rev. Lett. 97, 081303 (2006) 

[4] E. Babichev, V. Mukhanov, and A. Vikman, JHEP, 0802:101 (2008) 

[5] E. Babichev, V. Mukhanov, and A. Vikman, JHEP, 0609:061 (2006) 

[6] C. Eling and T. Jabobson, Class. Quant. Grav. 23, 5643 (2006) 

[7] D. Garfmkle, C. Eling, and T. Jacobson, Phys Rev. D76, 024003 (2007) 

[8] M. Choptuik, Phys. Rev. Lett. 70, 9 (1993) 

[9] J. Ziprick and G. Kunstattcr, Phys. Rev. D79, 101503 (2009) 

[10] R. Akhoury, D. Garfinklc, R. Saotome, and A. Vikman, in preparation. 

[11] G. Kunstatter, C. Leonard, R. Mann, and J. Ziprick, in preparation. 

A Maximal Slicing 

Here we present the details of the maximal slicing method. The metric is chosen to take 
the form 



which is the form for a spherically symmetric metric where the radial coordinate r is length 
in the radial direction in a constant time slice. Note that the usual area radius R is not one 
of the coordinates and is instead a function of the coordinates t and r. Thus the spatial 
metric ^y a b has components 



ds 2 



N 2 dt 2 + (dr + f3 r dt) 2 + R 2 (d6 2 + sm 2 9dip 2 ) 



(A.l) 



Irr = 1 



lee = R 2 
7 W = -R 2 sin 2 6> 



(A.2) 
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The extrinsic curvature, K ab is denned by 

K ab = -7a C V c n b (A.3) 

where n a is the unit normal to the surfaces of constant time t. However, due to spherical 
symmetry and maximal slicing, there is only one independent component of the extrinsic 
curvature. Specifically we have 

K\ = K*^ = -\K T r (A.4) 

Equation (A.3) is equivalent to 

d aij = -2NK tj + Dip, + Djf3i (A.5) 

where Di is the covariant derivative of the spatial metric 7$,-. The rr component of eq. 
(A.5) yields 

d r /3 r = NK\ (A.6) 

whose solution is 



P r = / NK\dr (A.7) 
Jo 

The 99 component of eq. (A.5) yields 

d t R = f3 r d r R + — RK r r (A.8) 

We now use the momentum constraint of the Einstein field equation to determine the 
extrinsic curvature. For maximal slicing (K = 0) this constraint is 

D a K ab = -K 7 6c n d T cd (A.9) 

Define the quantities P and S by 

P = ?i a V a (j), S = d r (f) (A.10) 

Then the quantity X is given by 

X = \{P 2 -S 2 ) (A.ll) 

while eq. (A.9) becomes 

R- 3 d r (R 3 K r r ) = -kC x PS (A. 12) 

The solution of eq.(A.12) is 

K r r = -kR~ 3 [ R 3 C X PS dr (A.13) 
Jo 

Note that there is also a Hamiltonian constraint associated with the Einstein field 
equation. In the case of maximal slicing, this constraint is 

(3) i? - K ab K ab = 2KT ab n a n b (A. 14) 
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where fi'R is the spatial scalar curvature. This equation yields 



d r d r R 



1 - (d r RY 



2R 



jR 



{K\f + 2k{C x P 2 - C) 



(A.15) 



We now determine the lapse N. It follows from the maximal slicing condition that 



D a D N = N 



K ab K ab +--T ab (n a n b + 1 ab ) 



(A.16) 



which yields 



8 r d r N +-(d r R)(d r N) 
Jrt 



N 



(K r r ) 2 + K(\C x {P 2 + S 2 )+C) 



(A.17) 



We now consider the evolution of the k-essence scalar field. From the definitions of P 
and 5 it follows that 



d t <j) = NP + (3 r S 
dtS = N(d r P + K r r S) + Pd r N + (3 r d r S 



(A.18) 
(A.19) 



The k-essence equation of motion, eq. (2.4) becomes after some straightforward but tedious 
algebra 



8 t P = p r d r P + Sd r N 
NC XX 



+ 
+ 



C x + P 2 C X x 

NC X 
C x + P 2 C XX 



{K\PS 2 + 2PSd r P - S 2 d r S) 
R- 2 d r (R 2 S) 



(A.20) 



The initial data are chosen as follows: P and are chosen freely, with S set equal to 
d r (f>. Then X is determined using eq. (A. 11) and then used to find C and Cx Then eq. 
(A. 13) is integrated to find K r r . Then eq. (A.15) is integrated for R using the fact that at 
the origin R = and d r R = 1. Then at each time step, the evolution proceeds as follows: 
first X is determined using eq. (A. 11) and then used to find C,Cx and Cxx- Then eq. 
(A. 13) is integrated to find K r r . And then eq. (A.17) is solved for the lapse N (using 
a tridiagonal method and the fact that d r N = at the origin and N — > 1 at infinity). 
Then eq. (A. 7) is integrated to find the shift j3 r . Finally, the quantities R, (ft, S and P are 
evolved to the next time step using eqs. (A. 8), (A.18), (A.19) and (A.20) respectively. The 
evolution is done using the iterated Crank-Nicholson method, and all spatial derivatives 
are found using standard centered differences. 
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